Geometric phases under the presence of a composite environment 
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We compute the geometric phase for a spin- 1/2 particle under the presence of a composite en- 
vironment, composed of an external bath (modeled by an infinite set of harmonic oscillators) and 
another spin- 1/2 particle. We consider both cases: an initial entanglement between the spin- 1/2 
particles and an initial product state in order to see if the initial entanglement has an enhancement 
effect on the geometric phase of one of the spins. We follow the nonunitary evolution of the reduced 
density matrix and evaluate the geometric phase for a single two-level system. We also show that 
the initial entanglement enhances the sturdiness of the geometric phase under the presence of an 
external composite environment. 
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I. INTRODUCTION 

It is widely known that a system can retain the infor- 
mation of its motion when it undergoes a cyclic evolution, 
in the form of a geometric phase (GP). This was first put 
forward by Pancharatnam in optics [l| and later studied 
explicitly by Berry in a general quantal system Q . Since 
then, great progress has been achieved in this field as it 
became clear that geometric phases had important con- 
sequences for quantum systems. The original notion of 
Berry phase has been extended to the case of non adi- 
abatic evolutions Q. As an important evolvement, the 
application of the geometric phase has been proposed in 
many fields, such as the geometric quantum computa- 
tion. Due to its global properties, the GP is propitious 
to construct fault tolerant quantum gates. The study of 
the GP was soon extended to open quantum systems. In 
this context, many authors have analyzed the correction 
to the GP under the influence of an external environ- 
ment using different approaches The interest on 
geometric phases has lately reached composite systems 
such as bipartite systems. In a previous article jl2l |. we 
have computed the GP for the bipartite system under the 
influence of a nonunitary evolution, induced by a bosonic 
or fermionic environment. There, we have shown that en- 
tanglement plays an important role in the robustness of 
the GP for open systems. In |13j . it was shown a partic- 



ular situation in which an initial separable state remains 
separable so that the GP of the system is always equal to 
the sum of the geometric phases of its subsystems under 
an unitary evolution. 

The presence of an environment can destroy all the 
traces of the quantumness of a system. All real world 
quantum systems interact with their surrounding envi- 
ronment to a greater or lesser extent. As the quantum 



system is in interaction with an environment defined as 
any degrees of freedom coupled to the system which can 
entangle its states, a degradation of pure states into mix- 
tures takes place. No matter how weak the coupling 
that prevents the system from being isolated, the evo- 
lution of an open quantum system is eventually plagued 
by nonunitary features like decoherence and dissipation. 
Decoherence, in particular, is a quantum effect whereby 
the system loses its ability to exhibit coherent behavior. 
Nowadays, decoherence stands as a serious obstacle in 
quantum information processing. 

In this context, there are many efforts to find a GP 
that "survives" the presence of the environment. What's 
more, the GP for a state under nonunitary evolution has 
not yet been directly measured. Recently, using a to- 
mographic approach, the GP was measured for a qubit 
coupled to a critical environment using a nuclear mag- 
netic resonance (NMR) quantum simulator [l^ . 

In this paper, we study the geometric phase (GP) of a 
qubit (spin 1) in a noisy 'composite" environment. That 
means, we start with a two-qubit-state (spin 1 and spin 
2) coupled to an external reservoir and then, we focus 
on only one qubit (by tracing out the other spin 2) in 
order to study the geometric phase (GP) of the resulting 
spin 1 (our subsystem). In the end, the model is dif- 
ferent to other preexisting models, such as for example, 
two-qubits in a noisy environment or/and a solely spin 
in a noisy environment. The reason for choosing this 
model is twofold. On one side, it is the natural sequel to 
other previous studies of the GP in open quantum sys- 
tems and has interesting features that differentiate this 
model from the others studied previously. On the other 
side, by taking into account new results about the ex- 
istence of an environmentally-induced GP for quantum 
states |l5|, we present a new insight on the effect that 
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the degree of entanglement has on the geometric phase 
of a particle of spin- 1/2 and, a possible way to measure it 
using NMR quantum simulators. Following this line, we 
shall start by describing the model and computing the 
geometric phase for a general environment in Section HIl 
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In Section Hill we shall focus on a bosonic environment 
writing explicitly the decoherence factors involved. We 
shall consider different situations for the initial state of 
the bipartite: either entangled (Sec. IIII Ap or a prod- 
uct state (Sec. IIII B[) . Further, in Section fill CI we shall 
study the differences that arise in the former computa- 
tions if one of the spin- 1/2 particles is not coupled to the 
external environment. Finally, we conclude our results in 
Section HVl 

II. DECOHERENCE INDUCED ON ONE 
TWO-LEVEL SYSTEM DUE TO A COMPOSITE 
ENVIRONMENT 

We shall consider a bipartite system (as in [12 )> that 
is to say, two interacting two-level systems, both coupled 
to an external reservoir. We shall study a model where 
the Hamiltonians that describe the complete evolution 
are given hy H = Hs + Hj + Hb where the system's and 



By F,, with i = 12,13,14,23,24,34, we arc implying 
the decoherence and dissipation factors that appear due 
to the presence of the environment in the reduced den- 
sity matrix. These factors are obtained from a noise and 
dissipation kernels similarly as in the familiar spin-boson 
treatment. The master equation for the reduced density 
matrix of a bipartite system coupled to a bosonic envi- 
ronment (Eq|3]) and the calculation of these terms can 
be read in the Appendix A of [l2j. However, we are not 
interested in the explicit form of these coefficients right 
now, since we want to focus on the geometric phase ac- 
quired by only one spin of the bipartite subsystem. For 
this, we need to know the dynamics of this only spin ruled 
by its own Hamiltonian and the presence of the compos- 
ite environment: the reservoir plus the other spin- 1/2 
particle of the subsystem. It is important to note that 
both spin- 1/2 particles interact by the coupling constant 
X defined in Hg and each of them interact as well with 
the external environment by the coupling constants Ai„ 
and X2n, respectively. We could further define the coef- 
ficients of the initial bipartite state so as to be able to 
interpret the results in terms of the known Bloch ball 
(we have not done so before because of the difficulty of 
writing Eq.Q). Therefore, we can assume that 

a = v/Aq cos(0o/2) 

(3 = -VA7sin(0o/2) 

C = V^sin((?o/2) 

6 = v/Arcos(6'o/2). (5) 



interaction Hamiltonians are respectively defined as 

Hs = -^<jI + ^tI+x Clival (1) 

N N 

Hi = CT]0^Ai„g„-l-cr^(8)^A2„g„. (2) 

n=l n=l 

The environment's Hamiltonian Hb will be left without 
definition until we are faced with a particular example. 

One approach to know the nonunitary features induced 
by the environment in the bipartite subsystem is by trac- 
ing out the degrees of freedom of the environment as 
suggested in Ref. [l^ . For a general arbitrary two-qubit 
initial state 

|«'(0)) = a|00) + I3\01) + C|10) + (5|11), (3) 
the reduced density matrix takes the form: 



(4) 



I 

The normalization of this state implies Ai = 1 — Aq. The 
concurrence of this state is C = 2-v/AoAi. Hence, we see 
that when Aq = Ai = 1/2 we have a maximally entangled 
state (MES). The angle Oq has here a geometric interpre- 
tation: coordinates X = (1 — 2Ao)sin(0o)j Y — and 
Z = (2Ao — 1) cos(0o) of the state vector in the Bloch 
ball are given, using the results obtained in [l3]- The 
radius of the Bloch ball is i? = |2Ao — 1| and all states 
lying on its surface have the same degree of purity (and 
come from two qubit states having the same amount of 
entanglement). The angle Oq is the one that makes the 
density matrix vector with the z axis. 



In order to know the dynamics of only one spin we 
need to trace out the degrees of freedom corresponding 
to the spin 2. After doing so, in this familiar notation, 
the reduced density matrix Pr{t) for the dynamics of spin 
1 is 



/3a*e^(2x+02)t^*^ |^|2 
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(Ao - 1 /2) cos 6^0 + 1/2 1/2 sin Ooe-^^^^Vit) 
1 /2 sin 6loe*"i*r(i)* (1/2 - Aq) cos 6*0 + 1 /2 



(6) 



r 



It is important to note that the T{t) factor includes a 
real decaying decohcrence factor as well as an imaginary 
dissipation term, both induced by the external parties 
to our subsystem spin 1. Therefore, we can write this 
coefficient as T{t) = r{t) exp^i-d^t)), where r{t) is the 
modulus of the complex number and its argument. 

The GP for a state under nonunitary evolution has 
been defined in as 



*9 = arg 5]v'£fe(0)efe(T)(vI'fc(0)|M'fc(r)) 



(7) 



where ek{t) are the eigenvalues and \'^k) the eigenstates 
of the reduced density matrix pr, and r is a time (r > 0) 
at which we study our system. Since we are working 
under the weak coupling limit, it is useful to consider a 
quasi cyclic path V : t e [0, r], with r — 27r/rii is the 
system's characteristic frequency). As we can see from 
the GP definition, in order to compute the phase gained 
by the system, it is crucial to know the system's dynamics 
at all times. This is why this is known as the kinematic 
approach in contrast to other approaches for GP in open 
quantum systems 0^ • The central result of Eq. ([7]) is 
to extract from the global phase gain acquired during the 



evolution, by a proper choice of the "parallel transport 
condition" , the purification independent part which can 
be termed a geometric phase because it is gauge invariant 
and reduces to the known results in the limit of an unitary 
evolution. 

The eigenvectors and eigenvalues of Eq.® are easily 
computed as 

e± = ±iy^cos2 6*0 + sin^ 9o + 4Ao(Ao - 1) cos^ 6*0 
1 

+ 2 



-init+i-i}{t) 



\v±) = 



r{t) sin ( 



r2 sin^ 6*0 + [(2e± - 1) + (1 - 2Ao) cos 6*012 

[(2£± - 1) + (1 - 2Ao)] 
r2 sin2 6*0 + [{2e± - 1) + (1 - 2Ao) cose'o]^ 



1) 



Unfamiliarly to all examples done before, this time, 
for an arbitrary initial state, we see that £+(0) 7^ 1 and 
e_(0) 7^ 0. Then, there is contribution to the GP coming 
from both eigenvalues. The only real factors in Eq.Q 
are a ~ ^e+(0)e+(T) and (5 = ^£_(0)e_(T). After 
some algebra, we obtain the geometric phase gained by 
spin 1 under this non-unitary evolution in a time t 



4>„ = arg( ar+(r)e*'^+e"^ + /3r_(T)e*^-e''^ 
I 



(8) 



where we have defined 



r-{T)e"f- = 

LP = 

We have also used 



K(0)K'+(r)), 

{v-m^-{r)), 
M 

(fii + —)cos^e+{t) dt, 
at 

(ill H )cos'^0-(t) dt. 

dt 



r{t) sin Oq 



(9) 
(10) 

(11) 
(12) 



COS0±(t) = 

y^r2 sin^ eo + [(2e± - 1) + (1 - 2Ao) cos Oof 
so as to be able to write the eigenvector in an easier way: 
\v±) = e-*"i*+'^(*)cos0±(i)|O) +sin0±(i)|l) (13) 



III. BOSONIC ENVIRONMENT 

In this Section, we investigate the effect of a bosonic en- 
vironment coupled to the composite system of two spin- 
1/2 particles, whose Hamiltonian is defined as 



Hb = 



N 

E 

=1 



(14) 



One assumption we shall make is that the spectral density 
of the bosonic environment J(w) is a reasonably smooth 
function of cj, and that is of the form up to some cutoff 
frequency A that may be large compared to fii and il2- 
The spectral density function can be written as J (to) = 
7o/4 a;"A"~^e~'^/'^, where dimensionless 701 ^ A^ and 
702 ^ A2 [ill • We can consider an ohmic spectral density 
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for such an environment, particularly one that goes as Abs|4>£]/;r 
J(w) ^ w. In that case, the r(<) factor that appears 
due to the tracing out of the degrees of freedom of the 
environment and the degrees of freedom of the spin 2, is 



A. Numerical results for an initial entangled state 



0.6 - >^ 



r{t) = e-2Toiog(i+A t )^(2Ao - l)cosnRt-ismnRt) 

(15) 

with Qji = {2x — 7oA). Wc have so far considered 701 = "■'^ ' \ 

7o2 — 7o for the sake of simplicity. \ 



0.2 - \^ 

0.0 r 'V 
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We shall start by considering the isolated spin, ini- 
tially set up in a bipartite state, but with no external 
interaction, that is to say 70 = and x = 0- Hence, in 
this uncoupled situation of the bipartite initial state and 
after tracing out spin number 2, there is a decoherence 
factor r{t) = (2Ao — 1) for spin 1. This means that the 
initial degree of entanglement will affect the dynamics 
through the decoherence factor. It is also important to 
note, that when Ao = 1/2, we have a maximal entangled 
state (MES), and the decoherence factor is F = 0, killing 
all coherences in Eq.® as expected. In Fig[T]we can see 
the behavior of the geometric phase (that wc call in 
the case of an initially entangled state) in the case spin 
1 is completely isolated but initially belonged to an en- 
tangled bipartite state. It is easy to see that it does not 
correspond to the geometric phase of a unitary single spin 
(i.e. = 7r(l — cos 6*0)) due to the initial entanglement. 
However, it can serve as a reference for studying either 
the influence of the environment or the entanglement of 
the bipartite state. One important feature to note is the 
monotonic behavior of the geometric phase in this case. 




FIG. 1: (Color online). Geometric phase 4>_b as a function 
of the concurrence C and the angle 9o for the uncoupled case 
X ~ and 70 ~ 0. The phase is measured in units of tt. 



FIG. 2: (Color online). Geometric phase of one spin initially 
set up in an entangled bipartite state for the isolated case 
70 = = X as a function of Oq. Red curve for C — 1, and 
purple curve for C = 0.06. The phase is measured in units of 

TV. 



If the degree of entanglement is Aq = 1/2, we have a 
MES and the concurrence is C = 1, for all values of Oq. 
In Fig|5] we present the geometric phase for two partic- 
ular initial entangled states of different concurrences for 
a better understanding of the behavior. In this figure, 
we present a solid red curve for maximum entanglement 
and a dashed purple one for a very low initial entangle- 
ment. We can note that for a spin initially set up in 
a MES the GP is 7r/2. On the other side, we can see 
that for an "almost" initial separable state (C ^ 0), the 
value of the geometric phase for only one spin depends 
strongly on the value of 9o (purple dashed line in Figl2|) . 
In any case, an initially separable (product) state shall 
be further discussed in a following Section. 

Let's recall that for a bipartite state, we can not longer 
use the Bloch sphere to seek a geometric representation 
of that state. In [l9| it has been shown that a geometric 
representation of a bipartite state can be obtained by 
using a Bloch ball and a SOS sphere. Therefore, the 
0-^ total phase gained by a state is a combination of not 
| only the dynamical and geometrical phase, but also 
0-4 the topological phase. Similarly to one qubit states, 
0.2 MES also gain a total phase of tt (or utt) under a cyclic 
0' evolution. However, this phase is of topological origin. 
It is already known that for a qubit the total phase 
gained is tt (or utt) and it is due to a combination of the 
dynamical phase and the GP. Now, we are obtaining like 
the "partial" phase of that entangled state, i.e. the GP 
of one spin- 1/2 particle after tracing out the "external" 
degree of freedom of the other spin, both initially set 
up in an isolated bipartite state. As a MES can not be 
represented in the Bloch Sphere, i.e. it is reduced to 
the central point of it, the eigenvalues of the reduced 
density matrix for Aq = 1/2 are degenerate. By using 
the definition of the GP for the degenerate case proposed 
also in Q, which in this case reduces mainly to the 
same formulation multiplied by the degeneracy of the 



5 



eigenvalues, it is easy to obtain the value of the GP for 
this case, namely 7r/2. Another case that can be easily 
analyzed is when do ~ 0. Hence, the initial state reduces 
to |*(0)) = VAo|00) + VI - Ao|ll), and it is easy to see 
that the cos 6±{t) = tt. Then, doing some algebra, we 
obtain that the $£; = tt for all Ao 1/2. On the other 
hand, if we have 6q ^ tt, which is equivalent to having 
l*(0)) 



the initial state |*(0)) = -^l - Ao|01) + V^|10>, the 
GP for a spin is zero for all Aq 7^ 1/2. These two initial 
states mentioned are the Werner states studied in (l^ . 
Thus, we see that the GP of a spin initially entangled in 
a bipartite state, yields a similar geometric phase than 
the one obtained for the whole bipartite isolated system. 




FIG. 3: (Color online). The evolution of the system can be 
illustrated by the path traversed in the Bloch sphere in the 
case that 70 = and x = 0. Red curve for C — 1, pink 
curve for C = 0.95, purple curve for C = 0.8, and black curve 
C = 0.43. Time is measured in units of Sli. 



As we are focusing on the geometric phase of one spin, 
we can find the path traversed in the Bloch sphere by 
studying the several components of the reduced density 
matrix p^- The three-dimensional coordinates in the 
Bloch sphere are x = + Pr-ni V = *(Pri2 ~ Pr2i) s-i^d 
z = — pr-i-2- By the use of the Eq.®, it is easy to 
generate the trajectories of FigI3] Thus, we can see a 
point for a MES, and bigger circles corresponding to the 
trajectories of the states with smaller degrees of initial 
entanglement, i.e smaller values of C, when x = and 
70 = 0. 

If we consider that there is no bosonic environment 
(70 — 0) but the two spin-1/2 particles are coupled 
through the x constant, then the decoherence factor is: 



TxW = (2Ao 



l)cos(2x0-«sin(2x<). 



phase as well. In all cases, except for Aq = 1/2, the degree 
of entanglement and coupling constant x act as a source 
of noise and dissipation (through its real and imaginary 
parts) for the system particle. In Fig|3]we show the be- 
havior of the geometric phase of the spin coupled to the 
other spin. One important feature that can be noted in 
FigUis that the monotonic behavior is broken. The main 
reason for this, is the contribution of both eigenenergies 
<p and (p in Eqs. (|TT|) and (|T2l) and the presence of a dis- 
sipation factor. Another interesting feature is that for a 
maximally entangled state (i.e. Aq = 1/2, C = 1) the ge- 
ometric phase is similar to that of the isolated case (40% 
bigger) . The presence of the valley is due to the intrinsic 
dynamics between the spin introduced by the coupling 
constant x- 



|*E|/7t0.5 




which explicitly shows that the coupling between both 
particles not only affects the dynamics through the time- 
dependent decoherence complex factor but the geometric 



FIG. 4: (Color online). Geometric phase <I>b as a function of 
the initial degree the entanglement through C and the angle 
60 for 70 = and x = O.lSli. The phase is normalized by tt. 



As a further step in our study, we include the interac- 
tion with the external reservoir. In FiglS]we can see the 
behavior of the GP as a function of the degree of entan- 
glement Ao and the initial angle ^o for 70 7^ and X 7^ 0- 
The presence of a stronger environment will contribute 
to the symmetry break and non-monotonic behavior, as 
can be seen in our case for a weak environment in Fig|S] 
Once again, we can note that the "open" GP gets a higher 
value for maximal initial concurrence and some particu- 
lar values of Oq. It is easy to note the non-monotonic 
behavior of the geometric phase, which exhibits a local 
maximum and minimum. This characteristic feature of a 
energy exchanging process between system and environ- 
ment (composite in this case) coincides with that shown 
in Ref.[2l| for a single qubit under a dissipative inter- 
action. As can be seen, the stability of phase can be 
significantly improved via a proper choice of the initial 
state determined by ^o which may be crucial for the effec- 
tiveness of quantum computation. Another feature that 
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FIG. 5: (Color online). Geometric phase $_b as a function 
of C and 60. The value of the geometric phase is measured 
in units of tt. Parameters used: x = 0.1f2i, 70 ~ 0.02 and 
A — 2OQ1 . The choice of the parameters used in for graphical 
reasons in order to compare with those plots of the preceding 
section. 



can be noted in Figl5]is that, as indicated in [l^, in some 
cases there is "saturation" value which is a characteris- 
tic value for a given configuration of parameters that do 
not change much in time. For example, an initial angle 
^0 = tt/S ~ 1, gives a maximum value of GP for a quasi- 
cyclic period of time. In order to study the importance 
of the initial state, we shall analyze the behavior of the 
geometric phase for different initial states in the weak 
coupling limit. 




5 10 15 20 25 

FIG. 6: (Color online). The geometric phase for the entangled 
initial state "1>_b as a function of time. The dashed pink line 
represents the initial state with Ao = 0.2; the dotted dashed 
purple line Ao = 0.1 and the black dotted dashed line is for 
Ao = 0.01 , all curves for 9o = n/5. Parameters used: A = 20, 
70 = 0.02 and x ~ 0.1. The geometric phase is normalized 

with TT. 



The geometric nature is the main feature of this phase 



in unitary closed systems. However, in the open quan- 
tum systems framework, this characteristic is not that 
evident. There have been many studies for one qubit 
systems but none of them has been conclusive about this 
topic. A nice way to see that the change of the GP is 
associated to the path traversed by the state of system is 
to plot the path traced by the state in the Bloch sphere. 
Therefore, Figs. [5] and [7] are presented to study the geo- 
metric nature of the GP in open quantum systems. Both 
figures can be related by noting that states with a higher 
value of initial concurrence, have a smaller rate of change 
of the GP in time (Fig. ^ and a smaller change in the 
path traversed (Fig. [7]), at least for small values of 9o- 
This is consistent with the geometric nature of the GP, 
which might be rephrased as it depends on the path tra- 
versed by the state of the system and not on the dynam- 
ics. The latter feature can be seen in the Bloch sphere, 
as initially bigger values of Ao, imply initially smaller ra- 
dius of the spiral trajectory and a smaller change in time 
of the GP undergoing an external influence. For exam- 
ple, if we compare the dotted dashed black curve with 
the dashed pink one in FiglHl we can see that the latter 
has a smaller rate of change in time than the former one. 
Likewise, in Fig 17] we can note that initially the black 
spiral radius is bigger and hence the change of the path 
and in the geometric phase is notable (during the first 
cycle, the black one 25% meanwhile the pink trajectory 
has changed only a 9%; the purple one changes a 16%). 
Therefore, we note that, not only is it important a proper 
choice of the initial value of 9o but the value of the initial 
concurrence plays a crucial role in order to have a robust 
GP undergoing an external influence. 

In this model we have several environmental parame- 
ters such as 7o, A and %. We can investigate the effect 
of the strength of the external bath, assuming a weak 
coupling limit, as we show the behavior of the phase as 
function of the initial concurrence and the coupling con- 
stant to the external bath. In FigUl we can see how this 
bosonic environment affects the different initial states of 
the system. It is easy to note that initially stronger corre- 
lations seem to contribute to the robustness of the GP in 
the presence of external couplings x and 70 . The value of 
X adds oscillations to the dynamics of the system. Once 
again, we can note that the stability of the phase can be 
significantly improved via a proper choice of the initial 
state determined in this case by the value of Ao and Bq. 
As for this figure, we can enhance our statement that 
initially stronger correlated states are prone to being less 
influenced by the presence of the environment and hence, 
the phase change is smaller for those states. In addi- 
tion, this figure serves for a future comparison among 
the equivalent behavior of the phase if the initial state is 
a product state as we will show in the next Section. 

It is noteworthy to mention the possibility of imple- 
menting this model in a NMR-quantum simulator, sim- 
ilar to the one used in Ref. [ij]. In this case, we have 
measured the GP-corrections over a qubit induced by an- 
other qubit emulating a critical bath. We have obtained 
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FIG. 7: (Color online). The evolution of the system can be 
illustrated by the path traversed in the Bloch sphere. Pink 
curve for C — 0.91, purple curve for C = 0.71, black curve for 
C = 0.43 . Parameters used: A = 20f7i, 70 = 0.02, x = O.lSli 
and ^0 = 7r/5. 




FIG. 8: (Color online). The geometric phase, measured in 
units of TT, as a function of the coupling constant to the 
bosonic bath 70, and the concurrence C. Parameters used: 
A = 2OQ1, X = 0, 6*0 = 7r/3 and r = 2tv/Qi. 



this correction by measuring the nonunitary evolution 
of the reduced density matrix of a spin- 1/2 coupled to 
an environment. The experiments were done using a 
NMR quantum simulator, where we have emulated quali- 
tatively the influence of a critical environment using just 
a simple onc-qubit model. By adding stochastic fields 
and further spins, we believe we can quantum-simulate 
more realistic environments and couplings to the system. 
With the same idea, it is possible to experimentally test 
the present model in which the entanglement parameter 
Aq plays an important role. In fact, we believe that this 
configuration will raise a new possibility of measuring 
the GP-corrections implementing a tomography of the 



reduced density matrix. However, this experimental test 
is out of the scope of the present paper. 



B. Numerical results for an initial product state 

In this section, we shall see what happens with the 
geometric phase (set as 3>p) of the spin- 1/2 particle if the 
initial state is not entangled, namely |^'(0)) = |4'i(0)) x 
1*2(0)), where 



|vi'i(o)) = v/T^|Oi) + v^|ii) 

1*2(0)) = v/T^|02) + \/^|l2>. 



(16) 



In this case, it is easy to associate the new parameters 
with the coefficients a, /3, and S in order to re-obtain 
the corresponding new dynamical equations: 



c 

S 



(17) 



An important consideration is to write the initial state 
of the bipartite system with equal parameters as in the 
entangled case, in order to be able to compare both cases. 
Therefore, we can associate p to cos^ (6*0/2) and write 



|«'i(0)) - cos(0o/2)|Oi) +sin(0o/2)|li 
1*2(0)) = 71^102) +V9|l2>- 



(18) 



In this way, we can describe the initial state of the system 
as a function of an angular and radial coordinates. As 
before, we can start by considering the case in which the 
bipartite system is isolated, that is to say 70 = X = 0. In 
such a case, we note that the geometric phase obtained 
is similar to that of a spin under unitary evolution 
<^>Pg = 27r(l - p), or equivalently $^ = 27rsin2 (6*0/2), 
with mod(27r), as shown in Fig 1^1 

The next steps in complexity is to consider that both 
spin-1/2 particles are coupled to each other through the 
X constant and/or to the external bath through the cou- 
pling constant 70. In such cases, the decoherence factor 
for the spin 1 (after tracing out spin 2) is 



^ p-27nlog(l+A''t^) 



cos(ri/ft) + i{2q — 1) sin(ri/}f) 



(19) 

In this case, there are also two eigenvalues, but only one 
contributes since e_(i = 0) = 0. So, the only contribu- 
tion comes from the eigenvalue e+. As the decoherence 
factor is complex once again, namely Tp{t) = r(t)e*''j'^*-', 
the computation of the geometric phase is different from 
the unitary one, namely 



$p = ar^ 



;|v/<W«(0)|<(r))e^^''(*)|, (20) 
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FIG. 9: (Color online). Geometric phase for the product 
initial state $p as a function of q (Ao) and ^o, for A — 20fii, 
X = 0, and 70 = 0. 



and 

\vl) = cosel{t)e-'^'*e'^'-^'^\0) +sm9l{t), (21) 

with 



0P/^^ ^ v/p(l~p)rp(t) 
+ v/p(l-p)rp(t)2+p2 




0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 10: (Color online). Geometric phase as a function of 
time for different initial states. Dashed lines (superposed) 
represent different values of q and 9o = tt/5 and the solid line 
is for q = 0.05 and 9o = tt/3. Parameters used: x ~ O.lili, 
A = 2OSI1, 70 ~ 0.02. The choice of the parameters used in 
for graphical reasons in order to compare with those plots of 
the preceding section. 

In order to analyze the importance of the initial state 
in the robustness of the GP in open quantum systems, 



we shall make the same analysis than in the preceding 
section. In FigHO] we present the GP as a function of 
the time for difFerent initial states in the limit of weak 
coupling with the environment. It is easy to see that the 
real crucial parameter is the initial value of 0o which sets 
up the initial state of |^'i(0)). The same behavior can be 
seen in Fig [TT] for the path traversed in the Bloch sphere. 
Let's recall that in the preceding section there was a 
hierarchy imposed by the value of the initial concurrence. 




1.0 



FIG. 11: (Color online). Trajectories in the Bloch sphere for 
the same initial states of Fig[71 Parameters used: x = O.lf^i, 
A = 20ni, 70 = 0.02. 




0.005 0.010 0.015 0.020 

FIG. 12: (Color online). The geometric phase, measured in 
units of TT, for different initial states as a function of the di- 
mensionless coupling constant to the bosonic bath 70. Pink 
lines for q = 0.4, 9 = n/3; and black lines for q — 0.01, 
9 — 7r/3. Dotted and dashed lines imply that x ~ O.lfii 
while solid lines imply x = 0. Parameters used: A = 20r2i 
and T — 27r/Sli. 



Finally, in FigH^l we present the GP for different 
initial states (the same ones considered before) as a 
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function of the coupling constant 70, for x = and 
X ^ 0. We can see that the presence of the internal 
dynamics introduced by x reduces considerably the 
value of the GP in the cases considered. 



C. The second spin is not coupled to the 
environment 

There is another situation that we can study which 
may result interesting. It is the case where the spin 2 is 
not coupled to the external environment, that means that 
the spin 1 is coupled to an external environment through 



This expression is very similar to that obtained for the 
case when the spin 2 is coupled to the environment in 
Section Fill Al However, if we look at the r'^(t) factor 
of Eq. p^ we see that in the former case there was 
dissipation introduced by the coupling of the spin 2 
and the environment, while in this case there is not. 
This is a somewhat reasonable result since if the spin 
2 is effectively coupled to the environment, when we 
ignore those degrees of freedom (by tracing them out), 
we expect to have a loss of energy. On the other side, 
just by tracing out the degrees of freedom of one spin, 
we do not expect to have a significant loss of energy. It 
is important to note that we have chosen to trace out 
the spin 2 but it is irrelevant whether we trace spin 1 or 
spin 2. The results are similar if we study the reduced 
density matrix for spin 2 by tracing out the degrees of 
freedom of spin 1. 



IV. CONCLUSIONS 

In this article, we have continued with the analysis of 
the geometric phase in the framework of a quantum open 
system, as a natural sequel to some previous works done 
in the field. In particular, we have analyzed the geomet- 
ric phase of a spin one-half coupled to a noisy composite 
environment. Initially we have started with a bipartite 
state coupled to an external bath and traced out the de- 
grees of freedom of one of the spins in order to focus on 
the GP of only one spin. The richness of the model lays 
in the fact that the initial bipartite state could be entan- 
gled or not, by introducing the initial entanglement as 
another free parameter of the model. We consider it as a 
major goal to understand the role that the degree of en- 



the Ai„ constant and to spin 2 through the x constant. 
This leads to a redefinition of the decoherence factors 
that appear in the reduced density matrix Eq.(l4]). For 
example, F12 = F34 = 1 and i^is = F14 = i^23 = ^24 = 
Fi where Fi = £-270 iog(i-i-A * ) jg -(.j^g decoherence factor 
obtain after tracing out the degrees of the environment, 
namely 

F^(t) = e-27oiog(i+A^t^)((2Ao - l)cos(2xi) -«sin(2xt)). 

After tracing out the second spin, we get the reduced 
density matrix for the dynamics of spin 1, for the case 
where the initial state is initially entangled. 



(23) 

I 

tanglement plays in the dynamics of the coupled system, 
as well as in the correction to the geometric phase. 

We have calculated the decoherence factor for different 
initial states of our system of interest and, consequently, 
wc have analyzed the geometric phase corrections from 
the nonunitary evolution. We have mainly considered 
two major groups: initial entangled bipartite states and 
initial product bipartite states, in order to study the 
role of the entanglement in the changes suffer by the 
geometric phase. In each case, we have computed the 
different decoherence factors and numerically obtained 
the open geometric phase for different physical situ- 
ations. With the help of the study of the geometric 
nature of the GP, we have shown that the entanglement 
enhances the sturdiness of the geometric phase under 
the presence of an external environment for small values 
of 6q . States with smaller values of initial entanglement 
seem to be less robust than others. This is related to 
the fact that the radius of the path traversed is bigger 
and its curvature is hence more affected by changes 
originated by the presence of the environment. Another 
important feature shown is that MES states still remain 
"special" or privileged as we have proved in Ref . [l2j due 
to the topological nature of the phase of these states. 
As shown, in order to implement a measurement of 
the noise-induced corrections to the GP, it is crucial 
the choice of initial state. However, we have shown 
so far that not only is it important a proper choice 
of the initial value of 6*0 but the value of the initial 
concurrence plays a crucial role in order to have a robust 
GP undergoing an external influence. In addition, it is 
important to state that by considering the case in which 
the temperature T is finite, the main properties of the 
geometric phase remain similar. However, it can be 
expected that the value of the maximum diminish with 



~Pr{t) = 



(Ao - l/2)cos(6'o) + 1/2 l/2sin(6'o)e-"^i*F^(t) 
I/2sin(6lo)e**^^i*F'^*(i) (1/2 - Ao) cos(e'o) + 1/2 
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increasing temperature. 

Furthermore, by comparison between the entangled 
and the product initial bipartite state, we have seen 
that having an initial entangled bipartite state yields 
a bigger geometric phase than in the ease of having 
an initial product bipartite state. In the particular 
case in which the central spin (system of interest) is 
completely isolated, it is easy to see that the GP of 
the initially entangled state does not correspond to the 
geometric phase of an unitarily evolving single spin, i.e. 

= 7r(l — cos 00 ) (due to the initial entanglement). 
However, for an initial product state, we have found 
that the GP, for the uncoupled bipartite system, is 
similar to that of a spin- 1/2 under an unitary evolu- 
tion. This suggests that the degree of entanglement 
works as an effective coupling between the two parti- 
cles in the system, even though one of them is traced out. 

Finally, we have noted that the stability of geometric 
phase with respect to decoherence and dissipation is 
crucial for effectiveness of holonomic quantum com- 
putation. It is evident that the stability of phase can 
be significantly improved via a proper choice of the 
initial state determined by the parameters of the initial 



state (whether Aq for the entangled state and for 
the product state). Thus, we consider the better choice 
an initial entangled bipartite state as opposed to a 
product initial state and with the feasible higher value 
of concurrence. 

All in all, we are reporting about a possible scenario 
where the phase measuring in an open system is feasible. 
We claim that the best choice in order to measure the 
open system GP must be take into account the degree of 
entanglement of the initial state. Our model admits the 
possibility of preparing a MES, which is the best option 
to intent a measure. This conclusion is a byproduct of 
having this type of model, in which one can prepare a two 
particle initial state, but using just one qubit to perform 
the experiment. 
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